Derealization of Flux Lines from Extended Defects by Bulk Randomness 
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We study the derealization by bulk randomness of a single flux line (FL) from an extended defect, 
such as a columnar pin or twin plane. In three dimensions, the FL is always bound to a planar 
defect, while there is an unpinning transition from a columnar pin. Transfer matrix simulations 
confirm this picture, and indicate that the divergence of the localization length from the columnar 
defect is governed by a liberation exponent v± — 1.3 ± 0.6, for which a "mean-field" estimate gives 
i>± ~ 0.78. The results, and their extensions, are compared to other theories. The effects may be 
observable in thin samples close to H a \. 
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It is now well known that fluctuations can drastically 
change the nature of the Abrikosov phase of type II su- 
perconductors. Thermal disorder can melt the vortex 
lattice, forming a flux liquid at both low and high vortex 
densities |Q . Quenched randomness also leads to new be- 
havior , which depends upon its degree of correlation. 
Point defects, such as vacancies or interstitials, encour- 
age line wandering, and may lead to the formation of a 
glassy phase with a non-zero critical current, called a vor- 
tex glass ||. In the extreme low density limit of a single 
flux line (FL), analytical J3J and numerical || work has 
shown that this is indeed the case, though the existence 
of the vortex glass away from H c i is still unclear. More 
recently, the experimental creation of "columnar," or lin- 
ear, pinning sites (f| has inspired theoretical treatment of 
the resulting "Bose glass" phase H. Twin planes, planar 
defects which occur naturally as a type of grain bound- 
ary, may lead to a different phase @] . 

In this letter, we study the behavior of an individual 
flux line in the presence of both point disorder and a sin- 
gle extended defect (e.g. a columnar pin; see Fig. |l|). The 
free energy of a d-dimensional FL of length L interacting 
with an n dimensional extended defect is 



+ V K (x(z),z)-AS{x ± (z))}, (1) 



where x(z) E 5R rf_1 is the displacement of the flux line, 
xj_(z) 6 is the component perpendicular to the 

defect, e is the stiffness, <5(x^) is a short-range potential 
concentrated near x^ = 0, and Vr(x, z) is a random 
potential representing the point disorder, with 

<V R (x, z)Vk(x', z')) = a 2 8{* - x')5(z - z'). (2) 

For n = 1 or n — 2, this model describes pinning by a 
columnar or a planar defect, respectively. 

In d < 3, the wanderings of a single flux line (FL) 
are controlled by the random potential due to impuri- 
ties. While in higher dimensions, a line freely fluctuating 
due to thermal effects is stable to weak randomness || , in 
this paper we consider mainly situations where the ran- 
domness is strong enough to be relevant. In such cases, 
the behavior at long length scales is dominated by a zero 
temperature fixed point However, even at T = 0, 
there can be a non-trivial transition between flux lines 
localized or free from the extended defect ||. The zero 
temperature phase transition also governs the singular 
behavior at finite temperatures. The relevance of disor- 
der, and the dominance of zero temperature fixed points 
simplifies certain aspects of the calculation. The parti- 
tion sum over all thermally fluctuating configurations is 
dominated by the optimal path which minimizes Eqjl], 
now regarded as an energy. 

The behavior of the FL in the absence of the ex- 
tended defect, A = 0, is well-understood from extensive 



theoretical and numerical work The energy per 

unit length is self- averaging and approaches a constant, 
E(L)/L — > (E(L))/L = —E . Fluctuations in energy and 
transverse extension also grow with length (albeit more 
slowly), and are described by non-trivial power-laws 



SE ~ AL U and Sx ~ BL< 



(3) 



where for short-range correlated disorder, the exponents 
uj and C depend only on d and obey uj = 2£ — 1 JIo| . £ 
is exactly 2/3 in d — 2, approximately 0.61 in d = 3, and 
gradually reduces to 1/2 in higher dimensions §. 

First, we obtain a lower critical dimension by consid- 
ering if the delocalized FL is stable to infinitesimal pin- 
ning, A. If the FL is confined by the defect, following 
Nattermann and Lipowsky jll| the localization length, 
£, is estimated by minimizing the FL energy 



E(L, i) = -AL/l d - n + AL/t 



EqL. 



(4) 



The first term is the attractive contribution from the de- 
fect, while the second term uses Eq.|| to account for the 
energy cost of confining the FL into L/i 1 ^ regions of 
length i 1 '^ and width I. Whether or not the pinning 
term is dominant at large distances determines its rele- 
vance. A weak potential is irrelevant for d > d/, where di 
is the lower critical dimension, defined by 



(di - n)( = I - uj. 



(5) 



Using the above estimates of (, di = 2 for columnar 
defects (n = 1), while 3 < di < 4 for planar defects 
(n = 2). A single flux line in three dimensions is thus 
always pinned by a planar defect. When weak pinning 
is relevant, minimizing Eqji] yields a localization length 
that diverges for small A as 



A " 0± , where v\ = 



1 — uj — (d — n)C 



(6) 



For the planar defect, v\ k. 3.6. These results are also 
obtained by power counting in Eq.|l|, taking into account 
the rescaling of temperature. 

For d > di, there is a transition between a delocalized 
phase for small A and a localized phase for large A. At 
such a transition, the localization length diverges as 



(7) 



where 8 = (A — A c )/A c is the reduced pinning strength. 
This divergence is accompanied by singular behavior in 
the energy, 



((E S (L,9))+E L)/L 



ll — Oi 



(8) 



which defines a "heat capacity" exponent a. For the 
borderline dimension di = 2 for n = 1, numerical simula- 
tions indicate a depinning transition as the defect energy 
is reduced B. 
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To construct a mean field (MF) theory of derealiza- 
tion, consider a FL very near the transition point on the 
localized side. Such a line makes large excursions away 
from the defect, forming "bubbles" of typical size z and 
I ~ z*> in the longitudinal and transverse directions (see 
Fig. |l|). For large z, the energy of a single bubble relative 
to a pinned segment is 



£(z) = (A - E )z - Az u 



(9) 



The first term is the energy cost (per unit length) of 
leaving the defect to wander in the bulk; the second term 
is a typical excess energy gain available (at scale z) due to 
a favorable arrangement of impurities. Minimizing this 
energy gives the scaling z ~ 0~"n and £ ~ 6~ v± , with 

v\\{MF)= 1/(1 -w), and v ± {MF) = Q/(l-u). (10) 

Assuming that the probability of encountering a favor- 
able bubble is independent of its length z, the total num- 
ber of favorable bubbles is proportional to L/z. The 
energy of the FL thus scales as 

E(L) oc-E L+-£(z) ~ -{E + d)L. (11) 

z 

The corresponding heat capacity exponent, a = 0, sat- 
isfies a modified hyperscaling form, 1 — a = (1 — 
appropriate to a zero temperature fixed point. 

Wc know of no obvious way of determining the upper 
critical dimension for the validity of the MF argument. 
In low dimensions, the argument may break down in a 
number of ways. The energy of a bubble configuration, in 
which interior returns are not allowed, may be different 
from Eq|| Intersections of the defect and the FL become 
probable when the sum of their fractal dimensions, n + 
1/C) is greater than d, yielding a critical dimension d u , 
which satisfies 



(du, - n)C = 1. 



(12) 



For columnar defects, 2 < d u < 3, while for planar de- 
fects, 3 < d u < 4. The above result is certainly a lower 
bound for the true upper critical dimension, but it is im- 
portant to realize that for the thermal depinning transi- 
tion it incorrectly gives d u — n = 2, rather than the exact 
result of d u — n = 4. Other factors, such as a renormal- 
ization of the linear term in Eq.^| by smaller scale bubbles 
or correlations between adjacent bubbles, may lead to a 
higher d u for the disorder-induced derealization as well. 
Taken at face value, the above result implies that the case 
of a columnar pin in three dimensions falls in the mean- 
field regime. Application of Eq.^Ojthen gives v± « 0.78, 
using the numerical values for £ and w in d = 3 Q . 

We now describe two approaches that attempt to go 
beyond the MF treatment. In the "necklace model" ap- 
plied by Lipowsky and Fisher (LF) to the borderline case 
of n = 1 and d — 2 , the FL partition function is de- 
composed into configurations with all possible sequences 



of pinned and unpinned (bubble) segments, and calcu- 
lated by Laplace transforms. For thermal derealization, 
since the partial partition functions of pinned and un- 
pinned segments are known, this can be done exactly. 
In the presence of randomness, the partition functions 
are not known exactly, and furthermore are themselves 
random quantities. LF propose using pre-averaged forms 
that depend only on the exponents C and u>. Although it 
is unclear that such pre-averaging faithfully reproduces 
the desired quenched quantities, LF's results agreed with 
both replica @] and numerical jl4j treatments of that 
case. A straightforward extension of their treatment to 
general n and d, gives 



v±(NL) 



i-J-n)C forO<(d-n)C<l 
(s4c^T forl<(d-n)C<2 (13) 
C for 2 < (d - n)C 



These results are substantially different from ours in both 
low and high dimensions. In particular, for the columnar 
pin in d = 3, the necklace model gives u±(NL) w 2.8. 
(Hwa has obtained the results of Eq.|l3] from a different 
approach pSfl.) 

A recent preprint Jl6| by Kolomeisky and Straley (KS) 
treats this problem for n = 1 by the renormalization 
group (RG), and concludes 



v±(KS)=< 



forQ<(d-l)C<l-w 
( d -i K -i+„ forl-^<(d-l)C<2-2c (14) 
j£- for2-2a;<(d-l)C 



1— o>— (d-l)C 



These results coincide with ours in the first regime, and 
at high dimensions. However, we remain unconvinced 
about the identification of the upper critical dimension, 
and the exponents in the intermediate regime []l9|| . For 
the case of the columnar pin in three dimensions, Eq.|l4| 
gives v±(KS) w 1.5. 

To differentiate between these theories, we examined 
the problem numerically for columnar and planar defects 
in three dimensions by a transfer-matrix method, which 
locates the optimal path exactly in a strip of finite width. 
To enhance performance, we chose the z direction along 
the diagonal of a cubic lattice, with random energies on 
the bonds. The energies of optimal paths terminating at 
position x at height z obey the recursion relation 

E(x,z+1)= min {E(x',z) +e(x,x',z)-A<5(x_L)}, (15) 

x— x'| = l 

where e(x, x',z) is a random energy for the bond con- 
necting x and x' at height z, and <5(xj_) is an appropri- 
ate lattice delta function indicating when the FL is on 
the defect. The above recursive computation is polyno- 
mial in the length L, allowing us to examine very long 
lines (L — 2.5 x 10 6 lattice constants with a transverse 
width W = 250 lattice units). The z-averaged end-point 
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displacement from the defect settles to a constant value, 
(xj_) z = £(W) 2 . In the localized phase, £(W) converges 
to a finite value as W — > oo, while in the delocalized 
phase, it grows as £{W) oc W. By examining £(W) for 
several pinning energies A and widths W, we find clear 
evidence for a depinning transition from a columnar pin, 
but none from a planar defect (see Fig. ^), consistent with 
our prediction for d/. Log-log fits to power-law forms 
yield v\ = 2.3 ± 0.1 for n = 2 and u x = 1.3 ± 0.6 for 
n = 1 (see Fig. ||). Distance from the critical region may 
cause additional systematic errors. The curvature in Fig. 
U suggests that this is indeed the case for the planar de- 
fect, which may explain the discrepancy with the result 
of Eq]| Given the result for n = 1, it is tempting to 
rule out Eq.|l3|. However, we must caution that related 
numerical simulations (l4[ in d = 2, see an effective expo- 
nent of fj_ = 1 (the mean-field estimate), before reaching 
the true asymptotic value of v± = 2. Resolution of this 
issue requires more extensive simulations. 

The actual situation in superconductors is complicated 
by several factors. Both defects and FLs appear at finite 
densities and with lengths limited by the sample thick- 
ness. First consider a single FL in a random set of colum- 
nar pins. The FL can be unpinned from a particular 
pin by two mechanisms. When the localization length, 
£, becomes comparable to the typical defect separation 

— 1/2 

£d = p d , the FL can freely wander between pins. (It 
may still be collectively localized by randomness in the 
distribution of columnar pins at a much larger length 
scale.) The second mechanism, originally described in 
the context of the bose glass 17 , involves a long-range 
hopping mechanism in which a FL wanders to a dis- 
tant pin. The energy cost of such a move (oc L 1 / 3 ) 
is compensated by the fluctuations in pinning energy 
(oc L 1 ! 2 j£ b with b = 3/(2C) - 2). A long enough sample, 
L 3> ihop £ 2 - 8 , will have sufficiently many favorable 
regions to allow such hops. The pinning effects described 
in this paper, are thus applicable only for L -C Lhop and 
£ <C id- A finite density of FLs introduces another length 
scale £fl = P FL ■ The dilute bose glass for £fl 3> is 
subject to the same constraints, and is unstable to the 
above hopping mechanism |l7j for L ^ Lhop- In the over 
dense limit of ^fl *C £d, the hopping mechanism is no 
longer relevant. A remnant of the unpinning transition 
may still be observed for £ <C £fl ■ Similar considerations 
hold for the case of grain boundary pinning. 

Columnar pins and grain boundaries are very effec- 
tive sources of FL pinning in superconductors. However, 
point defects in the bulk provide a competing mechanism 
for unpinning FLs from these extended defects. We pro- 
vide analytical and numerical arguments that the FL is 
always pinned to a planar defect, but that there is a tran- 
sition in the case of the columnar pin. The requirements 
for observation of the transition and the localized phase 
can be satisfied close to H c %, although the size of this 



region in the cuprite materials is extremely narrow. It 
is hoped that more sensitive probes will be developed to 
investigate these and other subtle effects in this regime. 
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FIG. 1. A flux line localized around a columnar pin in three 
dimensions (d = 3, n — 1). White portions of the FL indicate 
when it is on the defect. 
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FIG. 2. £(W) for various widths and defect energies for 

(a) the columnar pin and (b) the planar defect. The bond 
energies were uniformly distributed integers between and 
4095. The region of rapid change near A c « 375 indicates 
the transition for the columnar pin, while no such feature is 
present in the case of the planar defect. 

FIG. 3. Logarithmic plots of the localization length I ver- 
sus the reduced pinning strength for (a) the columnar pin and 

(b) the planar defect. The different curves in (a) account for 
the uncertain value of A c . 
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